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ABSTRACT 


V 

This  is  the  second  report  on  the  systematic  differences  be¬ 
tween  the  AGK3  and  SAOC  positions  and  proper  motions.  The  general 
theory  of  formulating  such  differences,  in  three  different  circum¬ 
stances,  is  developed  and  the  method  of  (hemi) spherical  harmonics 
applied.  Magnitude  and  color  terms  are  treated,  through  teems 
quadratic  in  m  and  c,  by  extending  the  variable  space  via  tensor 
product  spaces.  The  statistical  significance  of  an  expansion 
coefficient  is  ascertained  by  the  F  test.  Also  included  is  a 
discussion  of  the  time  dependence  of  these  differences  and  a 
resolution  of  most  of  the  poorly  matching  stars  (Table  2  of  Part  I) 
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I .  INTRODUCTION 

In  an  earlier  Report3*  we  presented  the  classical  method  of 

analyzing  the  positions  and  proper  motions  of  stars  in  common  to 

two  star  catalogs.  It  was  applied  to  the  Smithsonian  Astrophysical 

2 

Observatory  Star  Catalog  (the  SAOC  hereinafter)  and  the  third  of 

3 

the  Catalog  der  Astronomischen  Gesellschaft  (the  AGK3  hereinafter) . 
The  motivation  and  importance  of  such  work  was  given  there  too. 

This  Report  explores  the  method  of  using  a  complete  set  of  ortho- 
normal  basis  functions  to  perform  the  same  task.  Brosche  '  was 
the  first  to  propose  and  do  this  in  astronomy.  Moreoever,  follow¬ 
ing  the  public  presentation  of  the  first  part  of  reference  1®, 
we  received  two  suggestions  for  further  work.  Dr.  Clayton  \.  Smith 
(U.S.N.O.)  advised  us  to  look  into  the  time  dependence  of  the 
differences  between  the  AGK3  and  the  SAOC.  This  is  presented  in 
Appendix  II.  Dr.  John  M.  Sorvari  (LL-M.I.T.)  proposed  some 
astrophysical  and  observational  selection  effects  which  might 
explain  the  results  of  the  magnitude  and  color  effects  found 
earlier.  Our  discussions  and  ideas  on  these  points  are  in 
Appendix  III. 

The  next  section  describes  a  general  procedure  for  expanding 
any  appropriately  smooth  function  of  right  ascension  and  declina¬ 
tion  over  an  arbitary  "rectangular"  region  on  a  spherical  surface. 
(All  questions  of  convergence  and  smoothness  are  intentionally 
avoided.  The  theory  of  generalized  Fourier  expansions  is  well 
developed  and  we  can't  definitely  answer  convergence  questions 
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from  numerical  data.)  We  rapidly  specialize  to  hemispherical 
harmonics  for  the  northern  celestial  sphere  since  the  AGK3  >  nly 
extends  southward  to  6  ■  -2? 5.  The  incorporation  of  magnitude 
and  color  effects  occupies  our  attention  next.  All  of  this  is 
especially  easy  to  implement  computationally  when  the  parameter 
distribution  (i.e.,  a,  6,  m,  and  c)  is  both  uniform  and  dense. 

The  theoretical  basis  for  the  analysis,  in  two  different  ways, 
is  developed  next.  Sections  HID  and  IIIE  may  be  skipped,  there¬ 
fore,  without  loss. 

The  last  section  presents  the  results  for  the  positions  and 
proper  motions  with  and  without  the  magnitude  and  coloi.  terms. 
Appendix  I  is  a  full  FORTRAN  source  listing  of  our  P^fx)  generat¬ 
ing  subroutine.  It's  an  adaptation  of  the  subroutine  LEGPOL  from 
the  University  of  Rochester  Computing  Center  (#309.2.507).  The 
programming  for  the  majority  of  this  Report  was  performed  by 
Sharon  A.  Stansfield. 
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II.  EXPANSIONS  OVER  "RECTANGULAR"  REGIONS  ON  A  SPHERICAL  SURFACE 
By  "rectangular"  region  we  mean  any  region  of  a  spherical 
surface  that  can  be  described  as  consisting  of  the  simply  connected 

interior  to  a  *  a  ■  Oy  and  6  *  6^,  6  =  6^.  We  also  have 

0  <  a,  <  aT,  <  2tt,  -tt/2  <  5_  <  6T,  <  tt/2  and  we  don't  pay  attention 

—  L  —  U  —  —  —  U  — 

to  the  strict  inclusion  or  exclusion  at  any  boundary.  This 
definition  will  allow  us  to  handle  all  star  catalogs  of  interest. 

A.  In  General 

Define  t  the  usual  spherical  harmonic,  by 

Y~™(a*6)  *  N«™p,,™<sin5>exP(ima)  (la) 

run  nm  run 


where 


N 


nm 


f  2r+l  (n-m) 1 ]1/2 
l  4 tt  (n+m)  I  J 


(lb) 


and  P  is  the  associated  Legendre  function  of  order  n,  degree  m; 
nm 


Pn”<X>  2nnl 


/  i  \  <■)  in/  2  <n^in  n  n 

=  i^-a-X2)  r„{x^-l)n 


dx 


n+m 


,  m/2  a'Vfx) 
■ <-l)m(l-x2)  n 


dx 


m 


(lc) 


Here  Pn(x)  *  PnQ(x)  t*1°  or<*inary  Legendre  polynomial  of  order  n. 
The  set  n  "  0,1,2,...;  m  **  -n,  -n  +  l,...,n  forms  a 

complete  set  of  orthonormal  basis  functions  on  a  spherical  surface, 
viz. 


®  +n 

£  Yrnn(ct'6)  Ynm(a* '6' )  =  « (a-a • ) 6 (sinfi-sinfi • )  (2a) 

n=o  m=-n  ““  nm 


da 


Yr„Ja'6>Y£*m'  (Ot,6)COS6d6 
nm  n  m 


nn'^ram' 


(2b) 


where  the  asterisk  denotes  complex  conjugation,  6 (u)  is  the  Dirac 

delta  function,  and  6  is  the  Kronecker  delta  function. 

P3 

For  any  (se*  comment  on  page  1)  function  f(ct,6)  we  can  write 


f (a, 6 )  =  £ 


n=o 


+n 

E 

m=-n 


f  Y  (a ,  6 ) 
nm  nm 


where 


f 

nm 


da 


tt/2 

r  ( 

-tt/2 


f (a , 6 ) Y*  (a,6)cos<$d6 


(3a) 


(3b) 
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Now  define  (a,5)  by 
run 


*nm(a'6) 


[(“u“a 


1/2 


L1  (siridy-Bin^' 


N 


run 


run 


r2sin6- (sin6t+sin6y>  "| 

L  siniy-sinSj^  J 


(4) 


exp 


Then  it  follows  from  the  properties  of  the  that  for  any 

function  f(a,<5)  defined  on  the  "rectangle"  oeta^oty],  6e  [6-^,6^ 
(with  any  square  bracket  replaceable  by  a  parenthesis)  we  can  write 


*«■<«>  -  E  E  £nm 


n=o  m=-n 


(5a) 


with 


nm 


?  U  fU 

-  J  da  J  f  (a,6)y*in(ot,6)cos6d6 
6. 


(5b) 


B.  Hemispherical  Harmonics 

Because  the  AGK3  is  a  Northern  hemisphere  catalog,  we're 
primarily  interested  in  Northern  hemispherical  harmonics.  The 
Southern  extension  is  straightforward. 
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Define  H.(a,(5)  by 
nm 

Hnm(a,6)  "  /2  NnmPnm(2ain6-l)exp(ina)  (6) 

Then  the  forms  a  complete  set  of  orthonormal  basis  functions 

for  the  Northern  half  of  the  celestial  sphere.  We  note  in  passing 

that  only  the  combinations  considered  above,  Eqs.  (4,6),  have  the 

requisite  properties  of  being  both  complete  orthonormal  basis 

functions  in  their  domain  of  definition  and  require  no  ad  hoc 

presumptions  concerning  the  nature  of  f(a,6)  in  regions  outside 

7 

of  the  rectangle.  Hence,  the  proposal  of  Bien  et.  al.  [i.e., 
their  Eq.  (17)  wherein  a  Fourier  series  in  right  ascension 
decoupled  from  a  Legendre  polynomial  expansion  in  declination 
is  put  forward]  is  somewhat  curious.  Contrarily  (cf.  section 
IIC)  their  use  of  Hermite  polynomials  to  accommodate  the  magnitude 
term  appears  to  be  overkill. 

C.  Magnitude  and  Color  Effects 

The  functions  f(a,<$)  we  shall  deal  with  will  be  one  of 


Aacos<5 

=  (aAGK3~aSA0C 

)  cos6saoc 

A<5 

=s6AGK3"6SAOC 

(7) 

Ap 

s^(Aacos6)2  + 

2-)1/2 

(A«rj 

6 


r  2  2ii/2 

or  one  of  AyQ cos$,  Ayfi#  or  Ay  ■  (Ayacos6)  +  (Ay^)  similarly 

defined.  Not  only  will  these  quantities  depend  on  position , 
they  may  also  depend  on  the  brightness  and  temperature  of  the 
stars.  Thus  we  would  write  for  (say)  AS ,  AS  (a,  6  ,m  ,c)  .  To 
incorporate  this  in  Eq.  (6) ,  without  loss  of  generality,  we  write 


00 

y  r 

4 

0  s«o «- 

(mU-fl,L)  (CU"CL 

2*  -  (m  L+fliu)  ' 

•  p 

2c“ (CL+CU) ' 

mU^L 

s 

CU"CL  J 

r] 


1/2 


AS 


rs 

mn 


for  me[m^,niy]  and  celc^Cy]. 


(8) 
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III.  NUMERICAL  CONSIDERATIONS 


A.  In  General 

It  is  clear  that  in  any  numerical  implementation  of  the 
above  one  will  truncate  the  series  at  some  upper  limit,  say 
ne[0,N].  A  standard  result  of  analysis  is  that  the  best,  in  the 
least  squares  sense,  fit  of  some  function  f(x)  to  an  arbitrary 
orthonormal  basis  |pn|  on  xe[a,b]  to  order  N,  viz. 

N 

f(x)~£  fnPn  (x)  (9a) 

n=o 

occurs  when  the  expansion  coefficients  fn  are  given  by 
rh 

f n  = J  f(x)p*(x)dx.  (9b) 


To  be  explicit,  the  representation  (9b)  minimizes 


R n  =/  lf(x)  "  Z  fnpn(x)!2dx 

a  n— o 


(9c) 


Given  this,  it  would  appear  that  we're  left  with  a  numerical 
quadrature  problem  as  the  f  in  Eq.  (9b)  are  independent.  This 
isn't  quite  true.  As  long  as  N  isn't  infinite  there's  the  possi¬ 
bility  that  some  part  of  f  isn't  being  modeled.  This  is  an 
incompleteness  problem  and  will  be  reflected  by  the  fact  that 
is  still  larger  than  the  variance  due  to  accidental  (as 
opposed  tc  systematic)  effects.  The  independence  of  the 


8 


|fn|  does  assure  us,  though,  that  the  unmodeled  component 
of  f  isn’t  being  projected  elsewhere;  it  simply  will  be  left 
out. 

Having  disposed  of  that  point  we  are  left  with  a 
numerical  quadrature  problem.  When  the  data  densely  spans  the 
range  of  x  any  scheme  that  will  accurately  integrate  a  2N'th 
order  polynomial  will  do.  Note  though,  that  the  f  will  not 
necessarily  be  stable  as  a  function  of  the  order  of  the  quadrature 
scheme  because  as  we  go  to  higher  order  (which  also  means  higher 
resolution  in  x)  ws.  may  be  discovering  finer  details  in  f(x). 

If  the  data  isn't  uniformly  dense  or  the  value  of  N 
isn't  a  priori  decided,  then  the  problem  becomes  more  complex. 

In  the  first  case  the  sparseness  of  the  data  means  that  a 
numerical  implementation  of  Eq.  (9b)  may  destroy  the  ortho¬ 
gonality  property.  With  this  occurring  while  we  vary  the  value 
of  N  the  independence  of  the  |fn j  vanishes.  This  can  be  handled 
with  a  Gram- Schmidt  orthogonalization  procedure  (section  HID) 
as  N  varies  or  by  ordinary  least  squares  (section  IIIE)  with  N 
fixed.  The  reader  may  want  to  skip  these  two  tangential  sections 
on  a  first  reading.  Below  we  detail  the  numerical  procedures 
employed  here.  Almost  all  of  the  statistical  questions  relevant 
to  this  problem  one  might  pose  are  sloughed  off  here  because  we 
feel  that  nhey're  immaterial. 
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B.  Hemispherical  Harmonics 

The  values  of  Pnm(x)  were  computed  using  the  code  in 
Appendix  I.  The  values  obtained  were  checked  against  reference  8. 
Pnra(x)  was  computed  using  64  bit  arithmetic  and  then  truncated 
to  32  bits.  Gauss-Legendre  quadrature  schemes  using  16  and  32 
base  points  were  employed  to  perform  the  declination  integrals  in 
Eq.  (5b)  with  <$L  =0,  6y  =  90®.  We  judged  that  the  analytical 
simplicity  of  hemispherical  harmonics  outweighed  the  4.4%  loss 
engendered  |(2?5  x  360°)  [4tt (1806/tt)  2]  =  0.022  J.  This  still 
leaves  us  with  130,000  stars  nearly  uniformly  distributed  over 
the  Northern  half  of  the  celestial  sphere.  The  base  points  in 
each  case  are  the  zeroes  of  PL(x),  L  =  16  or  32.  Their  mean 
separation  is  5?3  +  2?3  and  2?7  +  1?2  and  if  PL(x^)  =  0  then 

=  sin-1[(l  +  x^)/2].  From  Eq.  (lc)  we  see  that  if  N  is  fixed 
at  15  (the  value  we've  used)  that  the  highest  order  "polynomial" 
we  need  to  handle  is  30.  The  normalization  integrals  for  the 
|Hnmi'  n  =  •  •  •  *15*  111  =  ~n»  -n  +  l..-.»n  were  performed 

'  -5 

as  a  check  and  never  differed  from  unity  by  more  than  5  x  10 
Newton-Cotes  type  formulas,  even  at  A5  =  0.01,  were  not  accurate. 
We  also  would  prefer  not  to  prejudice  the  results  by  having 
different  weights  at  different  positions  but  have  no  choice. 

The  closeness  of  the  mean  separations  in  declination, 
plus  our  just  mentioned  bias,  led  us  to  the  trapezoidal  rule  at 
Aa  =  5°  and  2?5  (as  appropriate)  to  perform  the  right  ascension 
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integral  (a^  =0,  ay  =  360®)  in  the  hemispherical  harmonic 
version  of  Eq.  (5b).  Again  the  normalization  integrals  were 
used  to  check  the  accuracy.  The  worst  deviation  from  unity  was 
0.014  at  5°  and  0.007  at  2?5.  The  possibility  of  using  a 
Filon  type  formula  was  also  investigated. 

The  details  of  forming  (say)  A5(a,6)  for  the  numerical 
integrations  were  as  follows:  Fix  L  =  16  or  32  and  A  a  =  5°  or 
2?5  correspondingly.  Compute  the  and  form  the.  midpoints  of 

each  bin  (more  complicated  splittings  have  no  better  theoretical 
basis  and  don't  move  the  boundaries  significantly).  We  now  have 
a  grid  over  the  Northern  half  of  the  celestial  sphere.  Differences 
for  each  star  are  placed  in  the  appropriate  elements  of  an  array 
and  the  average  value  of  each  difference  separately  computed  for 
each  bin.  These  are  the  numbers  used  for  (say)  A6(a,6).  In 
addition  to  computing  the  mean  we  also  calculated  the  variance 
in  each  bin  (with  respect  to  the  mean)  and  then  used  these  to 
provide  an  estimate  of  the  variance  of  an  individual  value  of 

A<$  • 

run 

C.  Magnitude  and  Color  Effects 

Formally  the  problem  is  as  above.  Practically  one 
must  insure  uniformity  and  density.  As  we  have  photographic 
magnitudes  for  all  of  the  stars  (from  the  AGK3) ,  as  long  as  we 
stay  within  the  bulk  of  the  magnitude  distribution  and  don't  use 
too  fine  an  areal  grid  there'll  be  no  difficulty.  From  Table  9 
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of  reference  1  we  choose,  therefore,  mL  -  B^O  and  niy  =  ll^O  and 
don't  look  for  more  than  a  quadratic  term  in  m  .  For  the  color 
indices  the  situation  is  not  as  good.  We  again  restrict  ourselves 
to  at  most  terras  in  c2  and  use  c^  =  0*1*0 ,  cy  =  270.  A  4'  point 
Gauss-Legendre  quadrature  scheme  was  used  in  each  case  and  to 
minimize  potential  difficulties  with  sparseness  the  16  x  5°  scheme 
was  used. 


D. 


Gram-Schmidt  Orthogonalization 
A  standard  result  of  linear  algebra  is  that  if  one  has 
a  set  of  K  linearly  independent  vectors  ju^  k  =  1,2,...,K  then 
there  exists  a  set  of  orthonormal  vectors  that  spans  the 

same  space.  The  |e^j  are  defined  through  the  intermediaries 


(10) 
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This  construction  is  known  as  the  Gram-Schmidt  orthogonalization 
process,  the  normalization  occurring  in  the  last  step. 

Consider  Eqs.  (9)  when  N  is  variable  and  the  values 
of  f(x)  are  not  dense  on  [a,b].  Then,  since  the  practical 
version  of  Eq.  (9a)  is 

f(x*>  -  £  fnpn(x,)  +*t.  I-  1,2 . L  (11) 

n=o 

for  L  samples  of  f  where  is  the  noise,  we  can't  simply  obtain 

f  via  the  extension  of  (9b) , 
n 

f„  -  *<x,)Pj(x,) 

Brosche4  when  dealing  with  this  situation  suggested  applying  the 
Gram-Schmidt  procedure  to  the  |pn(x£)}  •  In  particular  set 

En  =  IPn<xi) »Pn(x2) ' *  *  * 'pn(xL>3  and  define  {an}via 


9o  “  Eo 

an  -  E„  -  £ 1  JnA'  anm  "  (l2a) 

"  0  laj 

then,  with  akR  =  1  and  A  -  |anm},  p  -  {En}.  q  -  {<Jn} 

En  *  £  awA>  or  p  -  Aa  (12b) 

m— o 
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The  matrix  A  is  triangular  and  non-singular  (because  the  linearly 
independent  | £n  }  span  the  vector  space).  The  inverse  is  given 
by 

%  “  bnA  <13a> 

m=o 

where 


n-1 

bkk  *  *■’  brun  -  '  £  antbZm  for  ”  ‘  n  <13b) 

t=m 

and  AB  =  I,  the  N  dimensional  unit  matrix.  Equation  (11)  is 
transformed  into 


l  -  |  Fna„  *  5 

n=o 


(14a) 


where  f  *  |f(x^),  f  (Xj)  #  •  •  ♦  *f  (x^)  J  ,  E  is  the  new  noise  vector, 
and|Fn^  are  the  expansion  coefficient  of  f  in  the  basis.  As 
the  |  3Ln  |  are  mutually  orthogonal  the  least  squares  computation  of 
F  is  trivial. 


F  = 
n 


Vi 

IsJ' 


(14b) 


and  the  original  set  of  |fn|  are  given  by 


N 


:  =  y  F  b 
n  ^  m 


m=n 


mn 


(14c) 
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When  f  is  real  we  have 


where  f ( j ;k) m  f (a ^ ,  6 ^ ;m^)  etc. 
f  =  (-l)p  fp_qr  80  that  all  of  the  unknowns  aren't  independent. 
When  SpR  is  rewritten  solely  in  terms  of  independent  expansion 
coefficients,  the  gradient  of  SpR  computed,  and  then  set  equal  to 
zero  one  derives  the  normal  equations.  They  turn  out  to  be 


/ v  \P(3r  _  v  \iPSr» 
'pqrfmn  ypqrv£mn/ 

=  L  f (j;k)C/ra(j)Pn(k) 
j,k  tm  n 

uPqr  _  upqr\ 
pqrvfmn  ^pqr^mn/ 

-  £.  f(j!k)S£m(j)Pn0t) 
J 


(16  b) 


Wh-re  fpqr  ■  Xpqr  +  iypqr.Ypq(j)  =  Cpq(j)  +  iSpq‘j)  <X'y-C'S 

real)  and 


XPV  "  Cpq(3)Pr(k)Cim(j)pn(k) 

J 

“p“  "  Spq(j)Pr(k)Si.n(J>Pn(k) 

J  9  * 

=  ?v  Cpq(j)Pr<lt)Sto(i>pn<k> 
J  #  K 


(16c) 
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IV.  RESULTS 


In  order  to  decide  which  of  the  (say)  A  <5^  [cf.  Eq.  (18)] 
are  important  we  have  used  the  F-test.  Kendall  and  Stuart^ 
contains  a  readable  discussion  of  the  topic  of  discarding 
variables  including  the  fact  that  there  is  no  purely  logical 
way  to  do  so.  We  have  presumed  that  the  lower  index  expansion 
coefficients  are  a  priori  more  probably  significant  than  are 
the  higher  index  ones.  In  every  case  the  90%  level  of 
significance  has  been  used.  Moreover  because  the  systematic 
differences  are  real  the  expansion  coefficients  satisfy 

)  * 

Por itions 

Table  1  lists  those  expansion  coefficients  for 
AacosS,  A6,  and  Ap  which  were  significant  at  the  90%  level  when 
each  of  these  quantities  was  expanded  in  the  form  (5) . 

B.  Proper  Motions 

Table  2  is  similar  to  Table  1  but  for  Ayacos6,  Ay^, 

and  Ay. 


A6 


ns 

mn 
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TABLE  1 


EXPANSION  COEFFICIENTS  FOR  POSITIONS 


Aacos6  (0?001) 


n 


m 


Coefficient 


0 

0 

0 

0 

2 

0 

0 

0 

0 

0 


-0.808 

-1.47 

8.17 

-5.86 

3.71  -2 . 08i 
2.74 
-3.52 
-1.83 
6.15 
3.55 


A6(0?01) 


0 

16.8 

0 

12.9 

0 

-1.47 

1 

-2.29 

-6. 38i 

0 

12.7 

1 

-3.34 

+2 . 56i 

0 

2.57 

1 

-3.65 

+3 . 91i 

0 

10.5 

Ap  (O'. ’01) 


0 

0 

0 

0 


22.8 

0.713 

9.78 

-7.25 
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TABLE  2 


EXPANSION  COEFFICIENTS  FOR  PROPER  MOTIONS 


Ayacos6 


n  m 

0  0 

2  0 

3  2 

4  0 

5  2 

6  0 

7  0 

9  0 

10  0 


C0?001  cent) 


Coefficient 

-5.17 

20.9 

4.07  -8.22i 
18.7 

5.40  -10. 4i 
-13.3 
-8.10 
3.34 
8.89 


AUg  (01'01/cent) 


0 

1 

1 

2 

2 

3 

4 

5 

6 
9 

12 

13 


0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 


2.38 

-5.13 

3.72  +9.81i 

22.6 

5.54  -14. Oi 

13.2 
26.1 

17.2 
27.0 

28.2 
25.3 

-23.4 


Ap  (0" 01/cent) 


0 

4 


0 

0 


7.03 

11.6 


C.  Magnitude  and  Color  Effects 

Tables  3  and  4  list,  in  the  sane  format  as  Tables  1 
and  2,  the  coefficients  significant  at  the  90%  level  when  either 
magnitude  or  color  effects  are  separately  included  for  Aacoa6 
and  A6  as  in  (15a) . 
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TABLE  3 


EXPANSION  COEFFICIENTS  FOR  POSITIONS  WITH  MAGNITUDE  TERMS 

Aacosfi  (0? 001) 


n  m  r 

0  0  0 

110 
3  0  2 

2  2  0 

3  0  0 

3  2  0 

4  0  0 

4  2  0 

5  0  0 

5  10 

6  0  0 

6  0  1 

7  0  0 

7  10 

8  0  0 

9  0  0 


Coefficient 

-1.25 

1.51  +1. 17i 
-3.84 

6.96  +1.83i 
-8.83 

4.76  -3. 42i 
3.70 

-2.60  +1. 24i 
-4.40 

0.770  -5. 17i 
-1.10 
-4.27 
-1.93 

-2.10  +2 . 96i 
7.94 
1.74 


A6  (Ol'Ol) 


0  0  0 

10  1 
110 
2  0  0 

3  0  0 

3  10 

4  0  0 

4  10 

5  0  1 

5  10 

6  0  0 

6  10 

7  0  1 

8  0  0 

9  0  0 

11  0  0 

12  0  0 

13  1  0 

15  1  0 


22 

.0 

14, 

.6 

2, 

.36 

+1, 

.  45i 

14, 

.7 

-3, 

.82 

-2. 

.51 

-9, 

.071 

15, 

.6 

-5. 

,12 

+2. 

.  59i 

7. 

,29 

-5. 

,24 

+6. 

,  Oli 

12. 

,3 

4. 

,  15 

-6. 

84i 

-4. 

,67 

11. 

,6 

20. 

2 

-7. 

47 

6. 

83 

5. 

85 

-0. 

879i 

-3. 

29 

+0. 

918i 
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TABLE  4 


EXPANSION  COEFFICIENTS  FOR  POSITIONS  WITH  COLOR  TERMS 


n  m  r 

0  0  0 

10  0 
110 
2  0  0 

2  10 

3  0  0 

4  0  0 

5  0  0 

5  10 

6  0  0 

6  3  0 

7  0  0 

7  1  0 

8  0  0 

9  0  0 

10  0  0 

12  0  0 


Actcosfi  (OfoOl) 


Coefficient 

0.495 
3. 21 

4.26  +2 . 69i 
5.95 

-0.235  +3 . 98i 

4.12 
-1.78 
-6.47 

0.985  -3. 21i 
-6.39 

-2.20  +3 . 6°i 
-2.23 

-2.25  +3. 32i 
3.14 

6.27 
5.46 

3.12 


A6  (Oi’Ol) 


0  0  0 

10  0 
110 
2  0  0 

2  10 

2  2  0 

3  0  0 

3  0  1 

4  0  0 

4  10 

5  0  0 

5  0  1 

6  0  0 

7  0  0 

7  10 

8  0  0 

9  0  0 

10  0  0 

10  1  0 

13  0  0 

13  1  0 


10.2 

19.6 

4.34  +2. Hi 
14.0 

2.78  -5 . 55i 
-3.80  -9 . 21i 
2.47 
-3.23 
10.4 

-6.26  +2 . 84i 
20. 0 
6.25 
6. 38 
-4.08 

2.22  -3. 16i 
6.73 
11.3 

2  45 

-4 !  72  -1 . 05i 
-4.94 

5.13  -0. 858i 


APPENDIX  Is  COMPUTING  P  (x) 

nm 

As  mentioned  in  the  Introduction  we've  used  the  University 
of  Rochester  Computing  Center  subroutine  LEGPOL  to  compute  Pnm(x). 
We  had  to  modify  it  to  adequately  cover  the  n  =  m  case.  It's 
also  inaccurate  for  very  small  x  because  it  raises  x  to  the 
n  -  m  power  which  can  easily  result  in  underflows.  Several 
thousand  values  of  pnm(x)  were  generated  for  n  =  0,  1,  2,. ...15; 
m  =  -n,  -n  +  l,...n  and  compared  to  the  values  in  reference  8. 

The  largest  discrepancy  we  found  was  1  unit  in  the  last  decimal 

O 

plaee  (i.e.,  <5  parts  in  10  ). 

The  argument  list  for  LEGPOL  is  N,  M,  X,  ANS ,  AGAMA.  Here 

N=n,  M=m,  X  =  x,  ANS  =  P  (x) ,  and  AGAMA  is  a  vector  of  the 

nm 

natural  logarithm  of  the  gamma  function;  AGAMA (I)  =  £n[r(l+l)]. 
This  is  used  for  normalization  purposes.  The  full  FORTRAN  source 
code  follows. 
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r 


1 

2 
* 

4 

5 

6 

7 

8 
9 

IP 
11 
1 2 

13 

14 
Lb 
18 
17 
1 6 
iy 
2  'I 
21 
22 

23 

24 
2b 
26 
27 
26 

29 

30 

31 

32 

33 

34 
55 
3fa 
37 
3tt 
39 
4  0 

41 

42 

43 

44 

45 
4  b 
47 

46 

49 

50 

51 
5? 

53 

54 

55 

56 

57 
56 

59 

60 


SUBMJUTlNt  LE6P0L  L  N»n  »  y  »  A oJS  •  AS  API  A  ) 
IMPLICIT  KtAL*«  (A-H»<W> 

UintlMSl-jN  «bAMA(3l) 

I F  » Pi  •  Nr"  .  0  >  t»0  TO  6U 
IF  «i\l.NE,U  )00  TO  60 
63  ANbsi.00 
2U  KtluKi«l 

bU  IF  Lx.lM^.OHiO  TO  b2 
IF  .1)60  TO  62 

IF  l  m.S|  .1)60  TO  62 
60  10  b.3 
62  TE«t'l  =  A.ljO 
Ksl 
AsU 

L  1=2*14-1 
00  1  l=l*Llt2 

IF  »  I )  *6T  •  0  ) 60  TO  3 

2  Fulv  =  l.i)0 

1>U  I  L  1  fl 

3  MSA+1,1;;) 

FL)1V  =  A 

1U  I  EK(l=Tt;«H*FfMUM/FJlV 
1  hsK  +  1 

IF  LX. ti  l.  0)60  TO  5 
IF  Lim-1-1  )41«bb«40 
bb  S>U'*«=A 
60  tu  42 

41  Synsl.no 
60  TU  42 

4  U  ISlM-l-i 
LI=i/2 

23  F|\ioi  l=2  *'  -1 
FUl V=U 
FI  =  1 

CO  V  =  1,,.'0 

As A  *  *  1 
S  y  1*1  =  A 

00  2c  r*  =  1  »  L  I 
LOV  =  -L'»V/ 

PHt-=H*(Fl-l,UU  )  /  L  (Futw4  2.00)*FNUM*X**2) 
AsPKt  * 

SUI*'-6U'“i  +  M*LOV 
F 1=F 1-2.U0 
Fy  l  V=F  ■  >1  V  +  2  .  uU 
22  F0Oi'l  =  F  .  f)U 

42  lsPi/2 
Ksc*i 

Asl.U0-A*X 
F 1=1.00 

IF  Lrt-1  )  7  f%  /d,  79 
76  60  10  71 
7  7  60  16  70 
79  00  1 5  L  1  =  1  *  1 
7b  F I =F 1 #  A 

IF  L  (PI-*  )  .  (ML  .  U  )  6l  10  7i 
7  0  A(J6=  |  fc.3r|*SUi''i*F  1 
60  TU  20 

71  Alv6=ltM,'|*SUii*Fl*li6Llr<f  <  ) 

60  TO  2>i 
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52 


61 

62 

63 

b4 

65 

66 
67 
e»fi 
65 

70 

71 

72 

73 

74 

75 

76 

77 
70 
79 


b  10lrrsr,-M 

IF  (AJiFF  oEU*  0)  60  Tp  53 
lUcsiOiFF/2 

IFi  UUIFF-2*XQ2J  .Nfc.U)  GO  70 

XNUMSj  #qo 

XUtN=1.00 
00  30  1  =  1 « 10 IFF 
50  X;\|Ui1=Xn(JH*(  IlUFF -1+1 ) 
uo  31  r=l»102 

31  X0t-ivsxuEN*2.*l*(2.*N-?#*1+l) 

A’\)3S  I  Eki1*XNU|v/XL,E^ 

lF<rfUUaoiFF«4),.\(t.0>flRlS=-ANS 
60  10  £0 
52  MNb=U, 

60  10  <»u 
35  XNS!\I 

ANb=UEyp(AbAMA(2*N41,.AtoAMA(N+l>)*(l.D0/2.00,** 

tKU 
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APPENDIX  II:  THE  TIME  DEPENDENCE  OF  <AGK3-SAOC> 

The  table  gives  the  average  values ,  for  those  stars  in  both 
catalogs  with  matching  BD  numbers  and  with  the  indicated  Ap 
restrictions,  the  values  of  Ap  and  its  dispersion  (for  the  mean 
equator  and  equinox  of  1950.0)  for  the  epochs  of  place  1950.0, 
1975.0,  and  2000.0.  For  Ap  £  10"  the  means  of  the  components 
at  the  three  epochs  are  given  too.  Obviously  there  is  a 
relentless  secular  increase  in  all  of  these  quantities.  In  the 
Ap  £  2'. '5  bin  the  average  slope  is  l'.'5/cent  -  not  a  very  good 
state  of  affairs. 
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TABLE  A1 


THE  TIME 

DEPENDENCE  OF 

<AGK3-SAOC> 

Ap 

<  275 

AP 

<  570 

Epoch 

<Ap>  oA 

Ap 

(orooi) 

<Ap> 

(07 

o 

O  Q 

t— *  £> 

w  T3 

1950.0 

469  296 

473 

316 

1975.0 

897  515 

937 

535 

2000.0 

1206  592 

AP  <  10 ’.’0 

1423 

849 

Epoch 

<Aacos<5> 

<A6> 

(07001) 

<AP> 

aAP 

1950.0 

6 

42 

483 

419 

1975.0 

9 

54 

948 

649 

2000.0 

11 

66 

1448 

919 

C94-606 
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APPENDIX  Ills  ASTROPHYSICS  AND  THE  MAGNITUDE  AND  COLOR  EFFECTS 


Table  A2  repeats  a  small  part  of  the  information  contained 
in  Tables  10  and  11  of  reference  1.  The  questions  are  "Is  the 
trend  reversal  in  Aacos6  in  the  last  m  and  c  bins  meaningful?" 
and  "Is  the  discontinuity  in  AS  in  the  last  m  and  c  bins  meaning¬ 
ful?"  We  argue  in  posing  these  questions  that  A6  for 
b’Po  <  m  <  10?5  (the  ra  values  are  from  the  AGK3  and  presumably 

-  pg  -  pg 

homogeneous)  and  for  -O^S  <  c  £  1?5  (  c  =  m  -ra  both  apparent 

v 

magnitudes  from  the  SAOC)  are  constants  and  that  both 
9(Aacos6)/3m  and  3(Aacos6)/3c  are  negative.  The  constancy  of 
AS  is  beyond  statistical  reproach,  the  fact  that  3 (Aacos6 ) /3m<0 
nearly  so,  but  that  3 (Aacos6) /3c<0  questionable. 

Before  attempting  to  answer  these  queries  we  observe  that 
c  > lra  implies  that  the  star  is  red.  Red  stars  are  intrinsically 
faint.  From  Allen9  a  B-V  (=rav  -  mpg)  of  1?5  implies  Mv  =  9™8 
(Mv  =  absolute  visual  magnitude  =  mv  at  a  distance  of  lOpc  with 
no  absorption)  for  the  main  sequence.  B-V  =  1™6  means  My  =  11?8. 
These  stars,  therefore,  can  have  an  mv  ^  10m  if  and  only  if  they 
are  close,  at  about  lOpc.  Therefore  there  should  be  a  very  close 
correlation  between  the  c  >  l^S  and  ra  >  10m  results.  This  is 

pg 

exactly  what  is  seen.  Moreover  these  stars  must  be  distributed 
all  over  the  celestial  sphere  (because  they're  close)  and  because 
if  they  weren't  we'd  know  about  it  (e.g.  Gould's  Belt  for  the 
blue  O  and  B  stars).  This  means  that  they  can't  all  (or  even  a 
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sizeable  fraction)  to  be  from  one  of  the  source  catalogs  of  the 
SAOC  (except  for  the  AGK2) .  Hence,  we've  really  discovered  a 
systematic  difference  between  the  AGK2  and  AGIO.  We  now 
turn  to  potential  astrophysical  causes  or  explanations. 

Because  these  stars  are  close  they'll  have  larger  secular 
parallaxes  than  does  the  average  star.  But  even  if  the  mean 
epoch  difference  between  the  AGK3  and  the  SAOC  (read  AGK2)  is 
fifty  years  this  amounts  to  a  milliparsec.  It's  not  a  systematic 
difference  in  the  plate  reduction  schemes,  as  long  as  they  were 
the  same  (to  first  order)  because  this  is  a  catalog  to  catalog 
comparison.  Although  the  presence  of  a  color  and  magnitude  terms 
(^0'.'25|ra  -  9*1*121,  see  Eichhorn10)  in  the  original  version  of  the 

AGK2  is  known,  the  AGK2  plate  measurements  were  re-reduced  to 
provide  AGK3  proper  motions.  Similarly  we  rule  out  other  plate 
reduction  difficulties.  This  exhausts  the  astrophysical  causes 
that  immediately  come  to  mind.  Another  point  is  the  smallness 
of  the  numbers  in  Table  A2.  No  position  is  given  to  better  than 
O'. '01  and  very,  very  few  are  truly  known  this  accurately.  Hence, 
could  not  Table  A2  be  showing  us  truncation  effects?  No,  not 
in  both  right  ascension  and  declination  only  for  faint,  red  stars. 
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TABLE  A2 


THE  MAGNITUDE  AND  COLOR  EFFECTS* 


m  Range 

Aacos6 

A6 

( or ooi) 

Ap 

Number 

8™0  : 

9?0 

25 

38 

457 

25736 

9.0  : 

9.5 

4 

40 

430 

23247 

9.5  : 

10.0 

-10 

35 

444 

25576 

10.0  ; 

10.5 

-17 

40 

459 

20800 

10.5  : 

11.5 

-2 

57 

509 

15429 

All 

3 

40 

465 

126630 

c  Range 

-0™5  : 

075 

8 

31 

437 

19821 

0.5  : 

1.0 

-2 

25 

440 

24233 

1.0  : 

1.5 

-5 

39 

454 

21746 

1.5  : 

2.5 

12 

82 

499 

12069 

♦Stars  in  both  catalogs,  with  matching  BD  numbers,  north  of  the 
equator,  and  with  Ap  <  2'.'5. 
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APPENDIX  IV s  STARS  WITH  LARGE  POSITIONAL  DIFFERENCES 


Table  2  of  reference  1  lists  several  hundred  stars  for  which 
Ap  was  in  excess  of  6".  We  have  checked  each  of  these  stars 
individually  by  looking  then  up  in  the  SAOC  source  catalog.  For 
the  huge  majority  there  is  no  large  discrepancy  in  position. 
Rather  these  stars  are  members  of  a  binary  system,  each  star 
having  the  same  DM  number,  with  the  AGK3  containing  both 
components  while  the  SAOC  only  has  one.  As  we  sorted  by  in¬ 
creasing  right  ascension  (and  clearly  it  would've  mattered  how 
or  if  we  sorted)  within  common  DM  numbers  and  used  DM  number 
as  an  identifier,  nearly  half  of  all  such  occurrences  would 
result  in  a  misidentif ication.  The  majority  of  the  remaining 
large  Ap  values  apparently  belong  to  variable  stars.  A  few, 
probably  less  than  10,  represent  actual  catalog  errors  (e.  g. , 
the  SAOC  right  ascension  for  SAOC  92316  =  GC  3456  =  BD  +15  500 
is  wrong) .  These  are  still  being  sorted  out. 
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